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PREFACE 



This Investigation was prepared during the period 
February - June 1952 at the United States Naval Postgraduate 
School, Monterey, California. 

The subject was suggested to the author by previous 
experience in maintenance of small high-speed diesel 
engines, during T/h Ich burnt pistons had to be replaced 
on several occasions. 

The author wishes to extend his appreciation to 
Assistant Professor E. E. Drucker for his interest and 
advice throughout the development of the work. 
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SUKIMRY 



The objective of this work was to establish by 
annalytical means the temperature distribution in an 
internal combustion engine piston. For this purpose the 
piston was divided into a disk and a barrel, for which 
boundary conditions were expressed for the exterior 
surfaces and a common conical surface. The problem was solved 
in tv;o parts t first, a solution was found for the case where 
the gas temperature was assumed representable by an average 
temperature; this solution was then modified to account 
for a sinusoidal fluctuation of gas temperature. 

The resulting temperature distribution, within the 
limits of the approximations made, is: 
for the disk: 




for the barrel: 
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CHAPTER I 



INTRODUCTION 

Upon analyzing a piston of an Internal combustion 
engine with a viewpoint to understanding its behaviour 
as a heat dissipating member of the engine, it is readily 
realized that the piston absorbs a certain amount of heat 
through the upper surface of its crown. This heat is partly 
stored in the body of the piston until it reaches its steady 
state temperature; at the same time there is heat dissipation 
from the outside surface of the barrel and rings and also 
through its inside surfaces. After the steady state is 
reached, all heat absorbed must be dissipated through these 
surfaces. 

According to 0. L. Adams (1), five percent of the heat 
liberated by the combustion of the fuel must pass through 
the piston crown surface into the piston. Of this heat, ten 
percent is dissipated by the lower surface of the crown of 
the piston, as stated by J. L. Hepworth (5). 

Radiation from the gases into the piston only takes 
place through a very small portion of the time cycle and 
it amounts, according to B. Plnkel (9), to only ten percent 
of the total heat being absorbed by all the metalic parts- 
forming the container for the hot gases; thus it may be 
seen that the heat absorbed by the piston itself through 
radiation will be a part of this ten percent depending on 
the relative areas of the piston crown wl th respect to the 
area of the rest of the combustion chamber. Throughout this 
work it vl 11 be assumed that all heat absorption by the 
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piston takes place by convection. Also all beat dissipated 
by the piston vi 11 be assumed to do so by a process of 
convection. 

The temperature of the gases throughout the cycle 
varies widely; the temperature curves for an engine may 
be obtained either experimentally by direct measurement or 
annalytically from knowledge of the cycle under which the 
engine is operating. In this work the temperature of the 
gases will be assumed, for the sake of simplicity, to 
follow a sinusoidal variation about an average temperature. 

The main objective of this work will be to determine 
the temperature distribution throughout the piston. For 
this purpose a set of boundary conditions will be prescribed 
vihich resemble as closely as possible the actual conditions 
under which the piston operates. 
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CHAPTER II 



STATEIVENT OP THE PROBLEM 

1. Description of piston. 

A piston T/ylth the characteristics shovm in figure 1 
was selected for analysis. The crown will be a flat disk 
of uniform thickness 1. The barrel will be a cylindrical 
sleeve of constant thickness b. 

The boundary surface of disk and barrel will be 
surface a-a as shown in figure 1. This surface will have 
properties common to both disk and barrel and this will 
enable some of the unknown quantities entering the questions 
to be found. 

2. Temperature scales. 

Figure 2 shows a sketch of the temperature pattern 
that a particular point of the disk is expected to follow. 
This figure is given mainly to show graphically some of the 
values being used throughout the development. 

Figure 3 serves the same purpose as figure 2 with 
relation to the barrel. 

The temperature of the gases v/ill be assumed to be: 

where n will be the frequency of the cycle. 

The temperature Ta surrounding the outside surface of 
the barrel will be assumed constant. 

The temperatures for the disk will be measured from 
the level designated as in figure 2 and will be expressed 
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for any point of the disk as: 



The initial temperature v.'ill necessarily be a constant 
and will be designated as . 

The temperatures for points of the barrel 11 be 
measured in the j> scale of temperatures as shov;n in figure 
3, where: 



The initial temperature will be a constant. 

It is realized that the value of Ta will initially be 
equal to tg, but it is impossible with any degree of ease 
to consider into the problem its initial rise from tQ to 
Ta and thus it will be assumed constant and equal to T^* 

After a solution is reached, the results will be 
reverted to the standard level of measurement of temperatures 
in degrees P. which is designated as the t scale of 
temperatures . 

3. Boundary conditions. 

The following set of boundary conditions has been 
selected : 

(a). At the lower surface of the disk, the rate of 
heat rejection is ten percent of the rate of heat absorption 
at the upper surface. Since the rate of heat conduction is 
a function of the temperature gradient normal to the surface 
considered, and if we assume the coefficient of heat conduction 
to be constant throughout the piston, we may state: 
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Figure 2 
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There will be a difference in the areas of the upper 
and lower surfaces determined by the relative values of 
Rq and Ri. In this v;ork we will neglect this difference, 
but it is recognized that it could be included by multi- 
plying the value given as 10 by the ratio of Ri to Rq 
squared • 

(b). At the upper surface of the disk the rate of 
heat absorption from the gases will equal the rate of heat 
conduction from the surface into the disk, or: 



(c)« At the surface a-a the temperature of both the 
disk and the barrel will be equal. Since these temperatures 
will be measured from different levels, this relation may 
not be used directly without Introducing undesirable 
constants . 

Therefore it is necessary to realize that at this 
boundary the isotherms from disk and barrel must coincide 
and furthermore they must be continuous functions of the 
variables involved. This may be stated symbolically r 
along surface a-a: 




'Jr ■ ^3“ 



and : 



^ . ill 



' Jn.'' 
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(d). The lower surface of the barrel will be assumed 
insulated, or: 




(e). The Inner surface of the barrel will be assumed 
insulated, or: 



ii ) ^ 0 



(f). At the outer surface of the barrel the rate of 
heat rejection from the barrel into the surrounding medium 
will be equal to the rate of heat conduction from the 
barrel into its outer surface, or: 




(g). At time s zero, a point: 
at a temperature » or: 
at 

r-. a X--0 y- ^ 



z y 0, r ~ 0, 

B = 



will be 
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CHAPTER III 



SOLUTION OP PROBLEM 



1» General solution. 

In solving the problem the disk and barrel v?ill be 
treated separately and the solutions matched at the surface 

a-a. 

Both the disk and the barrel involve heat transfer in 
a cylinder; the disk may be considered as a flat solid 
cylinder, T^ile the barrel may be considered as a hollow 
cylinder. Under these conditions the best solution of 
Fourier’s Lav; of Conduction of Heat will be that expressed 
in cylindrical coordinates: 



Assuming the variables to be separable a solution as 
found in Appendix I 1st 



Prom figure 2 it may be seen that the expected curves 
for temperatures will be of the form: 



for a particular point on the disk. 

At the same time a sinusoidal variation will be applied, 
which proceeding along the z-axis will be damped and out of 
phase; therefore our solution would have a form such as: 




( 1 ) 
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(S) 
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• 

nl lelEnaTcf dsad evlovnl le-iiad add bnB delb add ddoS 

blloe iBil B EB btiPbienoo ed ^Bni delb edd jiabnll^o « 
wollod fl 8B be‘it>b fenoo ed vBm Ie*n;BO add eX t ^TebnlX^o 
“io aoldnXoF deed edd enoldibnoo esorid 'lebnU .TebnlX-^o 
beceeiqre darid ed XXlw JsaH 'lo noldcwbnoO lo wad E»ielnfo'5 

: eedBnlb<yooo XBol'ibnll'^o nl 

It 



I 



^ tk I ^ 

"* «- 1 i't vt 
6 



“Tt 



EB noldnios a eXciB'iBqes ed od eeldfiliBV odd anltaneeA 

,:ei I zibneqqA fil bnnol 

r - V. ^ ^ . '1 ^ 

^ ^ j^ 4 ctuD\ s ^ ^ ^ 

bscfo9cix0 d“Brfct n3ds 3cf nxoi*? 

:riol 9du lo 9cf lLt\7 ssnx/d'B'isqniact riol 



.ifelb *dd no dnloq ‘iBluoldnaq s no*i 
, btilXiqB ecf X£lw noldBl'iBv Xabloaunle fi eaild soibb ©dd dA 

dto bna beqeBb ed XIIb elxe-s add snoXfi gnlbeeooiq dolriw 
:«B rici/E nnol b evsd bX;^o« wolduXoe luo enoleiedd jesadq 



CX 




z 

T/here cz is the phase angle and e accounts for damping. 



Examining the previous equation it may he seen that 
the variables z and T may not be separated as assumed in 
the solution found by equation (2), and thus this equation 
will not yield an exact solution. 

A direct solution of the problem was attempted by the 
use of Duhamel’s Theorem (3), but it was found impossible 
to match the result obtained for the disk with that for 
the barrel at the surface a-a. 

Under these conditions it was decided to divide the 
problem into one where the gas temperature was assumed to 
be the average temperature Tjn> and then proceed to super- 
impose on it the effect of the sinusoidal variation in gas 
temperature . 

2. Solution for the disk. 

From equation (2): 







Since 9 must be finite throughout the body: 

B = 0 



and therefore; 









( 3 ) 



Prom boundary condition (a): 
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s 3 » 9 



(£) 



t^boX ^ri^ doo.f^rjondd sdlnlx ^d deun ^ 9onl3 

G = a 

:9iol9i9rid bna 

r o + ^^j>- 

oV ^ ^ fvi 3 3 s Q 



:(b) noIdXcnoo i^iebnood raona 






AhH'-' 



X 



Prom equation (3),: 

je _ f , f., , 1 , , , 

• r e ^ C h ■fJjfctraij J 'J 0 C^'>~/ 



therefore : 









0. I 



then: 






Q z ^ Q. I Cf> ^ \ 






y,.' U 7 -!- / \ 

^ ®3 

n>4^-0.( -» 



(4) 



or: 



From boundary condition (b): 




then: 

/^jLw i) ^ 

■ / 



Jll 

k 



The solution of this equation for b may only be done 
graphlcallyj it will yield a series of values for b vhich 
v;ill be designated as b^. 

Equation (4) now becomes: 

0 : 2 — s e 

V\ ^ / 




(5) 



Prom boundary condition (g): 

't'iO /l,-<5 

Substituting this in equation (5): 

G • &o 
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.<0 - C 



Therefore ; 



- -4<vVi,:j-rp , 

$ « 2_ ® j ■* j 



K »/ 



( 6 ) 



3. Solution for the barrel. 
Prom Equation (2): 



(j> z M e ^ t-y a ( /'^)J 

Prom boundary condition (d)t 
therefore ; 

i_ ~ ^ Lo 

then? 

f -. ? ^ ^ (3 - J? 

Prom boundary condition (e): 

li) = <? 



.^ 4 , 






or t 



A/ 



sf, 



-- F 



therefore : 

j> = 6le 









Prom boundary condition (f): 



Ip) =-^At ft ~ 



(7) 



( 8 ) 



(9) 
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Solving this equation graphically for the value of f, 
a set of values fp are found that satisfy this equation, 
and equation (9) becomes; 

i ■■ ‘ ^ a 



KC/p^}] 



(10) 



4. Conditions at the surface a-a. 
Recapitulating : 



^ ^ ^ >c[aZ )t 



^ 3 / 









kL 



and : 



Prom boundary condition (c) at points along the 



( 6 ) 






surface a-a; 










^3" ' 


E 


(11) 
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^ 

Oa'' 


Ja" 


(12) 



For equation (11) it is easily seen that for n*l 
and for a particular point along the surface a-a; 



KfO 



^ Z-. 



f 



ms., 



e-> 



then for n r 5 ; 






•s, 



VU 2/ 



this suggests that as a method of approach to the 
evaluation of a, g, and P we assume; 

g = bm 



(13) 



( 14 ) 
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Acceptance of this assumption involves the necessity. 



in order to maintain the equality of the above equation, 
that : 

Tm = Sp 

where : m :: p 

Also, in order that the effect of time on the tempera- 
ture be equal for the same point on the surface a-a as 
approached by either equation (6) or (10), that; 

fp = a 

Applying this condition to equation (13), we find: 







-L.]fi 






(IS) 



where ; 



m = p 

This equation is valid for points along the surface 
a-a, in. ere: J>/ . 



and should yield the values of Pp« 

Now, the replacement of equation (15) into equation 
(14) will yield an equation in terms of trigonometric and 
Bessel functions of r, valid for values of r ranging from 
Rl to Rq. Since this equation may not be simplified to any 
appreciable extent, it has been considered more advisable 
to replace the values of z and r of a particular point 
along the surface a-a, such that; 

z = ^/2 = 

r r (Ki + Ro)/2 =. Ra 
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q = III 
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( 
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: 9*19 rk ^s-B 
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then we find r 
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(16) 



where: m s p 

If now in equation (12) for n * 1 this same point in 
the surface a-a is investigated; 






Ji ) Vi i/? 1 



(17) 



Equations (16) and (17) should give the same value 
for Pp. There is no way of proving that this will he so 
except hy working out a particular problem* If in so doing 
it is found that the values of Pp satisfy simultaneously 
both equations, this will then mean that the assumption made 
that: g = t>iji 

and the derived expressions: 



- Sp 

and : f p I a 

are justified. 

In case that the values of Pp found by means of one 
of the equations mentioned do not satisfy the other 
equation, there is still the possibility of introducing 
this value in the second equation and determining a 
relationship betv.een Ri, b, Lq, and 1 such that both 
equations are satisfied and for which the assumption made 
is valid. This v;ould introduce a limitation to the physical 
proportions of the pistons for which this development is 
useful. The relationship just mentioned would take the formt 
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5. Partial solution using Tm» 

On the light of previous conclusions, the solution is 



now: 
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(f> ^£_f_ Vp e + 

1 m - I if? « f 



and since: 
and : 

we find, for the disk: 



& - t - 

.- 1 - r« 



t ‘ 1?*.^ +^51 ^ 

V X4 - , ^ i , 



;» 'oo / / i ^ 

e 



Oc*v J ^ vPaCj^ '’- J 



and for the barrel: 



t= T. +E'E.f’r f y.((,>)] 

>»»/ fiy ^ 



6. Superposition of sinusoidal variation. 

In order to superimpose the effect of the sinusoidal 
fluctuation of the gas temperature about a mean temperature 
Tjjj, the above equation (21) in a simplified form: 

t - T^<.s -f- i 



(19) 

(20) 

( 21 ) 

( 22 ) 
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we find that: 



after a quarter of a cycle, or at 




therefore : 






If this were put into equation (21) it would give an 
answer corresponding to the situation where this new 
temperature would have been imposed on the piston since 
time zero. 

To avoid this, instead oft. 



use the following expression given by M. Jakob and G. A. 
Eawkins (7): 



Tq Is the amplitude of the oscillation at the surface 
of the disk; 



The above expression v;as given for a' thick plate sub- 
jected to conditions similar to those of this problem. Even 
if in this case the disk may perhaps not be considered a 
thick plate, it is felt that this equation will yield a 
fairly good approximation if it is considered that it is 
common knowledge that the effects of fluctuating surface 
temperatures do not penetrate to any great extent beyond 
the surface exposed. In this particular case this effect 
is further minimized by the fact that the fluctuation of 
temperatures is very rapid. 






(23) 




(24) 



where : 
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Prom equations (23) and (24) It may be seen that what 
happens at the surface at time zero will have its effect 
at a depth z with a time lag of: 

, nir 



The addition of the term expressed in equation (24) 
will affect in equal form both the equations (21) and (22). 
Since for the barrel the values of z will, with the excep- 
tion of its topmost part, be far greater than those for 
the disk, this effect will only be perceptible in its 
upper section. 

Including equation (24) into equations (21) and (22) 
we have as a final result: 

For -the disk:- > 

• /<Tk"' 

. / /jzr ) 
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For the barrel: 

t; e 



fT^ 
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Vvis. f ff : 



7. Conclvisions. 

It is realized that the results obtained are not rigorou* 
but it is felt that they constitute a fairly good approxi- 
mation to the actual temperature distribution in the piston. 

Further work might be done on the subject by solving 
an actual problem in order to check the practicability of 
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the method used to determine Pp, By actual experiment on 
an engine the overall results of this work could be tested 
for accuracy and the percentage of error, if any could be 
found . 
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APPENDIX I 



GENERAL SOLUTION 



To find a general solution to Fourier’s Lav; of Heat 
Conduction for the case v;here the temperature should be 
analyzed in the unsteady state for a body best described 
in terms of cylindrical coordinates, we use; 



where ; 



Assume : 



Jr ' ^ J 

t~- 3, T ) 



( 1 ) 



then : 
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( 2 ) 



Since r, z, and T are Independent variables, we may 



say; 
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ct'T' 

Equation (3) may be rearranged as: 
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(3) 

(4) 

(5) 



which is Bessel’s equation of the first kind of order zero. 
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Then: 
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( 6 ) 



Equation 



which yields: 



(4) may be rearranged as: 

C <jr» i> ^ ^ 



Equation 



which yields: 



(5) may be rearranged as: 
S'f' 

0 = £ e 
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( 8 ) 



But , since : t - 9 ^ ^ ^ ^ 



we find using equations (6), 

Jf" r 

t - £ e ^ 



(7), and (8), that: 
JoUx) 



In this equation the constants a, b, A, B, C, D, and 
E must be found so as to satisfy the boundary conditions 
of the specific problem involved. 
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